Analytic Constructions of General n-Qubit Controlled Gates 
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In this Letter, we present two analytic expressions that most generally simulate n-qubit controlled- 
U gates with standard one-qubit gates and CNOT gates using exponential and polynomial complex- 
ity respectively. Explicit circuits and general expressions of decomposition are derived. The exact 
numbers of basic operations in these two schemes are given using gate counting technique. 

PACS numbers: 03.67.Lx 



I. Introduction- The study of quantum computers has 
been developing very rapidly over the past years. It pro- 
vides exponential speedup in factoring ll[ , or square-root 
speedup in unsorted database search [2j. In the circuit 
model of universal quantum computer the unitary 
operation that completes a computation task is a series 
of gates on a fixed number of qubits. Any unitarygate 
can be constructed from a set of universal gates [J 0]. 
Using the smallest number of basic gates to construct an 
arbitrary unitary transformation is very important, not 
only for using less executing time, but also for resulting 
less errors. 

Complexity of circuit is measured in terms of the num- 
ber of basic gates, namely the one-bit gate and the two- 
bit CNOT gate. For a general 2™ x 2™ unitary matrix 
U with 4™ degrees of freedom, 0(4™n 2 ) elementary op- 
erations are needed in principle [5(. Later on, efficient 
schemes implementing arbitrary quantum gates have re- 
duced the circuit complexity to 0(4 n ) 0, 0, 3. They 
are achieved by using the QR decomposition 6(, or the 
cos-sin decomposition 0] . General scheme for decompos- 
ing an arbitrary gate is given in Ref. [H using numerical 
method. For some quantum information task, such as 
initialization, a more efficient scheme with complexity 
0(2 n n 2 ) was proposed 

C n (U) gates are typical n-qubit fully controlled-?/ 
gates that apply a unitary U to the target qubit if and 
only if all the first n — 1 control qubits are 1. Cir- 
cuits for C 2 (U\ C 3 (U) and C 4 (U) gates have been con- 
structed [3, El E EH . But for the general case with 
n > 5, the explicit construction is absent. In this Let- 
ter, we present two different construction schemes for an 
arbitrary C n {U) gate, one uses an exponential and the 
other uses polynomial number of CNOT and one-qubit 
gates. The polynomial complexity scheme is good for 
large scale quantum computing. The exponential com- 
plexity scheme prevails for a circuit with a small qubit 
number. In particular, they are analytic. These results 



are very appealing in designing quantum computer pro- 
gramming language, because it not only saves computing 
time for its construction, but also avoids errors in numer- 
ical construction because of error accumulation. 

77. Exponential Construction Scheme- First we intro- 
duce some notation. For a generic n-qubit circuit, its 
qubits are numbered from the top from 1 to n. A k (V) 
stands for a controlled-!/ gate with k control qubits and 
one target qubit, so C n (U) gate is equally represented by 
A™-i(C/) w hose n—1 control qubits positioned at the top 
and the target qubit at the bottom. Order of operations 
in an expression as well as in circuits are performed from 
left to right. 

Previous investigations gave explicit networks of 
C n (U) gates for n = 2, 3, 4. In this Letter, we present 
a general analytic scheme implementing C n (U) gates for 
arbitrary values of n and any unitary operator U. Firstly, 
we define two kinds of quantum gate-array blocks, the A- 
block and the £>-block as shown in FigQ] 
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FIG. 1: A and B blocks in repair section. The left part is a 
A-block, the right one is a B-block. 

The ^4-block is indicated as A m , where m = 1, . . . , n — 
3. Its qubit nodes involve qubits m, m + 1, m + 2 and 
n. The 73-block is labeled as Bp where 1 < i < j < 
n. Its qubit nodes involve qubits i, j and n. First we 
suppose the explicit gate-array components of C n ^ 1 (U) 
network has been known, then we give a general analytic 
expression. Our strategy for C n (U) network is a two-step 
procedure: basic section constructing in the left part and 
repair section constructing in the right part of the circuit. 
Basic section is obtained by combining C n ~ 1 (U) network 
with a control input that is the (n — l)-th line without 
any performance. The basic section of C n (U) network 

is indicated as C"- 1 . C"- 1 contains 2 n ~ 2 - 2 CNOT 
gates, (2"~ 2 - 1) number of A 1 (U) and A 1 (U t ) gates, 
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where V = U. Repair section is yielded by placing 
A m and B* gate-array blocks in an alternating sequence 
with respective number of 2"~ 4 . 

A /3-bit Gray code [14j strings {g a }, where a — 
1, . . . , 2 /3 is a palindromelike ordering with special prop- 
erty that the adjacent bit strings differ only by a single 
bit. We define a function 7(0:, /3) to represent the numer- 
ical value of the position where g a and g a +i differ. In the 
repair section of C n (U), the index m of A m block is defi- 
nite as n — 3, the index j of 73* blocks is definite as n — 1, 
whereas index i varies complying with a (n— 4)-bit binary 

Gray code strings sequence. Denote C k as a network ob- 
tained from C k (U) gate combined with n — k extra qubits 
positioned between its last two qubits. Carrying out the 
recursion, the following results are obtained: 



C 5 (U) = C A A 2 B\A 2 Bl 

c 6 (u) = cT 5 a 3 bIa 3 bIa 3 b 2 a 3 bI 

C 7 {U) = C li A i BlA i BlA 4 BlA 4 BlA i BlA i BlA i BlA i B l 



C n (U) =C n - 1 A n - s B^ZtA n - s B2Zi...A n - s B^_ 1 . (1) 
So a generic C n (U) circuit where n > 5 can be expressed: 

J P+3 



C n (U) = (7 4 Y[ Jl A^ +1 Bf a ' 0) 



(2) 
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Let C|,(J7) denote 



Given a unitary operator V, there must exist one-qubit 
unitary operations D, E, F and real number a such 
that DEF = I and e la Da x Ea x F = V. a x and G are 
unitary one-qubit operations corresponding to matrices 

1 
e ia 

a A 1 ([/) gate where qubit k controls the qubit k'. We 
rewrite Eq. (2) in terms of CNOT and one-qubit gates 
after certain gate counting: 

n-4 2 13 

C"(U) = C* n n F^+U<T x )C^ +3 (a x )Gl +3 EiC^(<r x ) 

0=1 a = l 

C^ + \a x )G 0+3 E n C^+ 2 3 (a x )C^ +2 (a x )Gl +3 El 

y/3+3 



Cr£/\a x )Cl( a ' 9 \cT x )G p+3 E n CZ +i5 (<Tx)F n , 



(3) 



where 

C 4 = 



I>» C n {a x )E n Cl (a x ) C\ (g x )G\ El C 2 n (a x ) C\ (a x ) E n 
Cl(<J x ) C\ [a x ) EIG\C\ (a x ) C\ (a x )G 3 E n C 2 Z (a x ) 
Cl(a x )GlEtCi{a x )C^{a x )E n C 3 {a x )GiG 2 G 3 F n . 

(4) 

In Eqs. (EE]), C%{a x ) are CNOT gates and the D, E, 
F, G and their hermitian conjugate are the one-qubit 
gates and their subscripts represent the positions. The 
/3-bit binary Gray code strings reflected in a sequence of 
BZ^^ can be chosen freely in an arbitrary cyclic /3-qubit 
Gray code sequence and the C n {U) circuit for (n > 3) 
are self-inverse. 



We can prove this exponential simulation fulfills the 
action of C n {U) faithfully. After a carefully accounting 
of merges of CNOT gates and single-qubit, we find this 
exponential simulation scheme for a C n (U) gate finally 
utilizes 2" — 2 CNOT gates and 2™ one-qubit gates. 

777. Polynomial Construction Scheme - The above 
scheme is advantageous for small values of n, but it be- 
comes inefficient for a large value of n, for instance n > 8 
because of its exponential complexity. Here we propose 
a C n (U) circuit using 0(n 2 ) basic CNOT and one-qubit 
gates. 

We know for n > 3, C n (U) can be simulated by a 
network with its own inverse, where V 2 = U in Fig. H 
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FIG. 2: A quantum circuit for implementing C n (U) g; 
with unitary V meeting V 2 = U. 



;ates 



Given the explicit construction of arbitrary C™ 1 (U) 
gate is known, the key procedure is to simulate two 



'{(J x ) gates. For n > 4 and m\ € 1, . 



1 — 2, a 

f\ n ~ 2 (cr x ) gate can be partitioned into two A mi (cr x ) gates 
and two A 1 ™ 2 (a x ) gates, where mi + 7712 = n — 1 as shown 
in Fig. [3] 
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FIG. 3: A quantum circuit for implementing a A n 2 (cr x ) gate 
with A mi (da,) and A 1 ™ 2 (er^) gates where mi + 7712 = n — 1. 

So the problem is reduced to how to construct A™ 1 (a x ) 
and A m2 (a x ) gates. If we assign mi = [n/2], ni2 
11 — [n/2] — 1, A'" 1 ((7a;) and A m2 (cr x ) can be decomposed 
into several Toffoli gates. It is worthy noting that these 
decompositions are only applicable to A mi (cr x ) for n > 6 
and to A m2 (a x ) for n > 7. So for n > 6, we investigate 
the most regular arrangement of Toffoli gates implement- 
ing A mi (a x ) in Fig. g] Let T c h denote a Toffoli gate with 
control qubits a and 6, the target qubit c. Consider the 
characteristic features of such network: A" 11 (cr x ) network 
consists of 4 [n/2] — 8 Toffoli gates and the indices of Tof- 
foli gates is symmetric around io = [n/2] — 1 and periodic 
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FIG. 4: A quantum circuit for a A 1 ™ 1 (a x ) gate consisting of 
Toffoli gates arranged most regularly. 



with period d = 2[n/2] — 4, then use a formula to describe 
the regulation of Toffoli gates 



A™ 1 (a x ) 



n r- 



4[n/2]-8 2+/(n) 

n l+(l-«« (i0 )(l-5 i , i0+<i )(n-[n/2]+/(n)) 
-n-[n/2]+2+/(n) ' 

1=1 

(5) 

f arctan(tan(f i - §)) - [§] + 1|. It 



1 when 1 < i < d 



where f(n) = 

denotes the deviation of i to 
or to ?j — 1 + d when d + 1 < i < 2d. The absolute 
value function expresses symmetric property, the arctan 
function fixes periodic regulations, and the 5 functions 
correspond to certain singular points at i = [n/2] — 1 and 
3[n/2] — 5 referred in above formalism. 

Similarly, for n > 7, /\ m2 (<J x ) gates can be simulated 



by the network in Fig. [5l We find that there are An 
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FIG. 5: A quantum circuit for a A m2 (cr2;) gate consisting of 
Toffoli gates arranged most regularly. 

4 [n/2] — 12 Toffoli gates in A m2 (a x ) network. Inspecting 
the mathematic property of indices, it is periodic with 
period d' = 2n — 2 [n/2] — 6, symmetric around jo — 
n — [n/2] — 2 and jo + d' = in — 3[n/2] — 8, we obtain the 
following formula 

4„-4[„/2]-12 ™+(l-^.j )(l-^, J - 0+d /)<2[«/2]-2n + 3 + g( „)) 
A 2 (<Tx) = [Q T n-l + ( S l-5 3 l )(l-5 3 2 „_ 2[ „ /2 j_ 5 )(2[»/2]-2n+5+ 9 (n»' ( 6 ) 

where g(n) = |f + £ arctan(tan(f j- f )) -n+ [§ ] + 2|. 
Then we propose a cascade decomposition of C n (U)(n > 
7) gate by a recursive method shown in Fig. [6l where 
unitary Vi is defined by = U. 
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FIG. 6: Cascade structure of C n {U) gates where unitary Vi satisfies V/ 2 = [/. 



Suppose r p^ 1 - 2 ' - k U i s a /\k 1 (a x ) whose first k — 1 
qubits control the last one qubit. Then above decompo- 
sition can be described in formula 

G n (U) = C™' r ^V' 1 )T^' 2 '-" _2} C™ _1 (V 1 +)T^' 2, '-"~ 2} C^ _1 (V'i) 
= CwT-e) n T<if -■ fc - 2} C^ 1 (vt_ fc + 1 )T{i. 2 '- " fc - 2 >c£- V„_ fc + 1 ), 

where j^b^-' 5 - 2 } = A™ 1 (<j a ) A™ 2 (cr x ) A™ 1 (a x ) A™ 2 
(<r x ). After tedious calculation, we find out a C n (U) 
gate totally requires 2 CNOT gates, 8n 2 - 72n +174 



Toffoli gates and 2n — 3 two-qubit controlled gates. 
Then the problem is reduced to simulating Toffoli gate 
with basic CNOT and one-qubit gates. Using well- 
known congruent modulo phase shift (CMPS) meth- 

ods M for Toffoli gates, it can be expressed as if = 
R c C b c {a x )R c C a c {a x )RlC b c {a x )Rl where R = R y (n/4). 
The CMPS scheme only requires 7 basic operations which 
is much less than 14 basic operations in the usual simu- 
lation scheme. C & (V n -o) part is congruent to the circuit 
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for C 6 (U) and we have proven that the Toffoli gates la- 
beled as 4, 6, 8, 10, 15, 20, 25, 30 as shown in Fig. 

for C 6 (V n -6) part, and all the Toffoli gates other than 

C 6 (V n -e) in Eq. (7) can be replaced by the modulo phase 
shift of Tollofi gates. 
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FIG. 7: The explicit structure of C 6 (U) (congruent to 
C 6 (V n -6)) in terms of Toffoli and two-qubit controlled gates. 



Taking account of the merges of CNOT gates and one- 
qubit gates, we obtain the total number of basic opera- 
tions in C n (U) construction are 24n 2 -212n + 540 CNOT 
gates and 32n 2 — 288n + 739 one bit gates ultimately. 



IV. Conclusions- In conclusion, we have given two an- 
alytic schemes for constructing a C n (U) gate for arbi- 
trary value of n and any unitary U operator, one with ex- 
ponential complexity and the other with polynomial com- 
plexity. General expression for decomposition of C n (U) 
gats with basic one-qubit gates and CNOT gates has been 
derived explicitly. We have compared the exact num- 
bers of basic operations required in these two methods 
for n = 1 — 20. It shows that the exponential construc- 
tion is advantageous for the value of n = 1 — 8, whereas 
the polynomial simulation is efficient for larger values of 
n > 8. 



Given a unitary operator V n ~k+i, L, P, R and S 
are one-qubit unitary gates such that e lb La x Pa x Q — 

V n -k+i, LPQ = I and S = * J 6 

represent which qubit they are performed on. Now we ob- 
tain the C n (U) in terms of CNOT and one-qubit gates: 



their subscripts 
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C n (U) = C6(V„_ 6 ) n Wk^Wk^st^QiC^^PlctH^Ll 
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